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A SPECIAL CASE OF THE LAPLACE CO-EFFICIENTS 4\°. 
By Pror. AsapH Hatt, Washington, D. C. 


FEBRUARY, 1887. No. I. 


(1.) When we consider two planets moving around the sun in the same plane, 
and take this plane for the plane of reference, the potential, or perturbative func- 
tion will be, in the notation of the J/canique Céleste, 
m'r cos (7! — v) 


r, rv’, v, v’ are the polar co-ordinates of the planets referred to the centre of the 
sun, and m’ is the mass of the disturbing planet. If the eccentricities of the or- 
bits are small, the radii vectores will not differ much from the semi-major axes of 
the orbits. Let these be a and a’, and suppose a < a’, a=a/a’; also let v’ — » 
=. Then we have to develop into a periodical series the function, 


I 
(1 —24acosé+ 


Omit the subscript s of Laplace, and assume 


= 16" + cos + cos 20 + 6” cos 


+ +... 


In the Mécanique Céleste the values of the 6” co-efficients are given in series 
which proceed according to powers of a. For the principal planets of our solar 
system this method is generally convenient, since the greatest value of @ in this 
system is 0.723, which occurs in the theory of Venus and the Earth. But for 
values of « greater than this the series converge slowly, and its seems better to 
determine these co-efficients by a different process. In the second volume of his 
great work on elliptic integrals, Legendre shows how the 6" co-efficients may be 
reduced to these integrals, and their numerical values computed by means of his 
tables. In this case, however, the iritegrals can be reduced to forms so conve- 
nient that it seems worth while to proceed directly, and without reference to tables. 
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If we integrate the expression cos 74 cos j# d#, between the limits o and =z, 


we have, 


f cos 74 cos 70 dt = (¢+7)4+ cos(¢ —/) 0 
Oo 


sin (2 sin (7 —j) 


re) 


Since 7 and 7 are whole numbers, the integral becomes zero, except for ¢ = 7, when 
its value is $x. If ¢ and / are both zero, the value of the integral is z. Denot- 
ing, therefore, the left-hand side of the equation for the series by (4), then mul- 


tiplying the equation by : cos 7 d#, and integrating from o to z, the terms are 


all zero except 


e 


xX cos 74 cos 76 di, 
oO 


We have, therefore, for the general value of the co-efficient, 


=x cos 70 di. (1) 


Thus for 70, and 7=1, we have 


 — 2 a 
— 2acos@ + a’) 


2 cos dd 


- 


Let us consider three successive co-efficients of the series, 


=? cos (i — 1) 
oO 
cos it db, 
Oo 


a 2 e 
git =? cos (i + 1) 0d0; 
Oo 
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we can easily find a linear relation between these co-efficients. Integrating by 
parts the expression 

cos — 2acos 4 + a’). dd, 
we have 


frcos it vy (1 — 2acos 4 + a’). dd 


(1 —2ucos a) 


f sin @ sin 20 dé 
y (1 — 2acos + a?) 


The first term is zero at the limits o and =, and since 
sin # sin = 4 [cos 1) — cos(é + 1) 9], 
the result is 


f cos (1 — 2acos 4 + a?) di = — 


ce) 


Treating the same integral in a different manner, we have 


py (1 — 2a cos + a?) dd 


cos 74 di cos 4 cos 74 


cos 7 do cos + 1)@d0 
(1 — 2acos 0 + a’) 


af cos (¢— 1)4d0 


By comparison of the two values we have 


2 
This equation shows that if we know the co-efficients 6'~" and 4", we can 


compute from their values the next co-efficient 6°+". Hence all the co-efficients 
may be computed from 4" and 6'”. We can find also from the definite integrals 
the derivatives of the 4 co-efficients with respect to 4, which are required in the 
development of the disturbing forces, and these derivatives may be made to de- 


if 
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pend on the 4" co-efficients; but these formula are well known and appear to 
be convenient for computation. 

(2.) In order that we may reduce the expressions for 4 and 6" to forms 
suitable for numerical computation, we have to transform the integral containing 
— 2acos # + as adivisor. Since cos = 2 cos? — 1, we have 


(1 — 2acos + = + 2a + — 4acos? 


(1 


If, therefore, we put 


(i + a)?’ 1+ 

2d. 
r=cos}#, — 


the value of 4°’ becomes 


ax 


This elliptic integral may be transformed according to the method given by 


- Lagrange.* For this purpose put 


Then Br?) = — 27) (1 — 
and ax x dx 

— #7) (1 — — 
Equation (3) gives — 2x? — 


and differentiating, we have first, 


x dx dy 
— 14+ — 227° 


Again, solving (3) for +”, we find 
14+ FP =y[1+2(F—2)¥ + 


so that dy dy 


* Oeuvres de Lagrange, Tome 11. p. 253. 
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Assume the trinomial under the radical to be of the form (1 + a”) (1 + ,3)*), 
and we shall have 
a= —2+ 2, (1—), 
—2—2y(1 —#); 
and the trinomial is resolved into the real factors, 
1+ 
The quantities a and 3 are of the same sign since their product is #*. The half 
sum of these quantities is 4* — 2, which is negative, and therefore a and / are 
negative. 
Let 


from which we have 


(1 —F). 
Thus we have 
ax dy 
Now we have 


and, if in equation (3) we change our unit, or put 


i+« 


2 (=a). 


our integral becomes 


ax =ittf dy 
J vy —4*) (1 — 2 — 


This last integral is of the same form as the first one, and it is plain that we 
may continue such transformations as far as we please. Hence if we put as be- 


fore 
fit 
— ay 


“= 


etc. ; 
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we shall have 


VL — +?) (1 — 


17H 1+% 1+4 f dy 


From the values of a, a, 4, . . . &, we see that 4, can be diminished at 
pleasure, and, at last, a,?y? may be neglected. In this case we have 


From equation (3) we see that y is zero when + = 0, and when + = 1; and 
that y has a maximum for an intermediate value of +. Hence in the first step of 
the transformation we must double the integral for the limits, y=oto y= 1; 
and the same for each succeeding step; so that in the final value we shall have 
the factor 2". The result will be | 


= “a, 1 Oy I+ 4, I+ 4, 
or +4). . .(1-+4,). (4) 


The value of 4" may be tound by means of the same substitutions, which 
give 


I 
aS [i — 24) (1 — 
The preceding reductions give 
27 — 1 = —y[it+ 2(F—2) + 
and making the same change of the unit as before, we have 


20? — 1 = 1) + fa — — — 
Putting, therefore, for dx / ,/[(1 — (1 — its value found before, we have 


(247 — 1)dx _ « _ (a —1)dy I+a 


2 


2 


i 
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It is evident that this method of transformation may also be continued, and 
at last we shall have the forms 


(2x? — 1) dr 


J — 4%) (1 — 


I I 
I+a@i+a 
oO oO 
a +4)... (1+4,)) 
a.a,.(1 I 1 
| 


The first term approaches zero as its limit and may be neglected. The integrals 
of the second set vanish at the limits zero and unity. We have left only the last 
terms within the brackets, and since as before the value of the last integral is 
2", we have finally, 


Equations (4) and (5) can be computed easily, and then the remaining co- 
efficients may be found from equation (2), when @ is not very small. 

We may notice here a relation between the moduli of the successive trans- 
formations which was pointed out by Lagrange, and which is similar to that em- 
ployed by Gauss in his celebrated memoir on the secular perturbations of a 
planet: Determinatio Attractionis, etc. 1818. 


From the values of a, 4, a, etc. we have 


and hence is greater than g, greater than etc. 


q 

= <=. cfc.; 
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We havealso 


etc., etc. 
If we put pt+q=m 
etc., etc. ; 
so that p=km+n), 
A=Hm+m), Hm — mM); 
we shall have mM =m+n, N= 2,/(mn); 
My = + My = (mM); 
etc., etc. 


In the series of quantities, therefore, 172, . Mh, My, . ., the 
corresponding terms are the arithmetical and geometrical means of twice the 
terms that precede them. 

(3.) For the principal planets of our solar system there does not appear to 
be much need of departing from the methods commonly given for computing the 
4" co-efficients ; but in the Saturnian system we find a case where these methods 
are not convenient. The two satellites of Saturn, Titan and Hyperion, the largest 
and the smallest of this system, move in the same plane, and from the values of 
their mean distances from the centre of the planet we have the ratio 


0,8250863. 


Hence the usual series for computing the 6" co-efficients will converge slowly, 
and the preceding method by definite integrals is much easier. I will compute 
the values of some of these co-efficients for this value of 4. For this purpose put 


a=sing, 4 %=sing,, ete. 


The expressions for @, 4, . . . give 
I —cos¢ 
= — tan? Le, 
I + cos ¢ 
I —cos¢ 
— — tan 
I cos gy 
I — cos Co 
ad, = = }¢,, etc.; 
I + COS 
and hence I + a, = sec’ hy, 


1 
2 

I a, = sec? 1y,, 
L 
2 


I + = etc. 
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The computation of 5°” and 4 by formule (4) and (5) may be arranged as fol- 
lows : — 
a= sin ¢g = 0.8250863, 


log sin g = 9.9164993760. 


Calculation of 6. 


¢ =55° 35/ 50%07 log log 
dg = 27 47 55.035 tan 9.7219831665 sec 0.0532568941 


¢,=16 8 16.349 sin 9.4439663330 
4¢,= 8 4 8.174 tan 9.1515781796 sec 0.0043209827 
I g 5.8228 sin 8.3031563592 


O 34 32.9114 tan 8.0021702287 sec 0.0000219319 
¢3= O O 20.8337. sin 6.0043404574 


4g, = 104168 tan 5.7033093 _0.0000000005 
0.0575998092 
2 log sec = 0.115 1996184 
log 2 = 0.3010299957 
log J = 0.4162296141 
= 2.6075 31805 


Calculation of 


log @ = 9.9164993700 = 0.825086 3000 
log = 9.4439063330 
log } = 9.698970004 3 2nd term = 0.1 14666277 3 
2nd term = 9.0594357133 
. log a, = 8.3031563592 3rd term = 0.0011522908 
log = 9.698970004 3 
3rd term = 7.0615620768 4th term = 0,00000005 8 2 
log 4, = 6.0043404574 0.940904025 3 
log = 9.6989700043 log = 9.9735457423 
4th term = 2.7648725 385 log $ = 9.698970004 3 
log 4, = 1.4066186 log 6” = 0.4162296141 
log } = 9.6989700 log 6 = 0.088745 3607 


5th term = 3 8704511 — 20 == 1,226719759 
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Having found 6 and 4", we can compute the remaining co-efficients by 
formula (2), Thus 


I+ @ 
b 5 . 3.6 


etc., etc. 


The preceding calculations would be much shorter if made with seven-figure 
logarithms, and the accuracy would be sufficient for practical purposes ; but the 
work has been done with ten-figure logarithms in order to show better the ra- 
pidity of the convergence. The value of 4'" computed by the common series, 
including twenty six terms and a, is 

= 2.6075 285 3. 

The following table contains the values of the first six co-efficients com- 
puted by means of the definite integrals and formula (2); and also the values of 
the same co-efficients computed from the formule of Leverrier.* 


Values of a= 0.8250863. 


Co efficient. by Value by series. 
2.6075 318 2.6072929 
1,2267198 1.2265599 
0.7967742 
0.5624428 0.5622774 
0.4129396 0.41269C6 


0.3102719 0.3 100983 


The values computed by Leverrier’s formula include twelve terms, but they 
are too small, and it would be necessary to continue the series to double this 
number of terms if we wish to obtain values accurate in the sixth decimal place. 
But if series are to be used, they should be transformed so as to proceed accord- 
ing to powers of I — a’, as was pointed out by Euler in 1750. 

(4.) The history of these co-efficients is interesting. Since it is necessary 


* Annales de [ Observatoire de Paris, Tome IU. Addition premiére, p. 5. 
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to consider the reciprocal of the distance between two bodies when their mutual 
attractions are to be determined, and the best practicable method of finding these 
attractions is generally the expansion of this reciprocal into a periodical series, 
mathematicians were immediately confronted with the problem of finding the co- 
efficients of the sines and cosines of this series. D’Alembert (/ean /e Rond) saw 
in 1743 that all the 4’ co-efficients could be found as soon as the first two were 
known, a relation which is expressed in equation (2). Euler and Clairaut made 
many improvements in the actual numerical determination of these co-efficients. 
The elegant method of resolving the general expression of the reciprocal, 
(1 — 2acos @ +- a?) ~*, into two binomial factors by means of imaginary quanti- 
ties, which is employed by Laplace, is due to Lagrange, It is probable that the 
transformation of the elliptic integral which has been used above was suggested to 
Lagrange by the work of Landen. It has since been employed by several writers, 
and among others by Professor Woodhouse.* 

In looking at the great amount of labor that has been done on the impor- 
tant problem of developing the perturbative function, one cannot fail to notice the 
variety of the symbols that have been introduced. This variety has become so 
great as to cause confusion, and in order to be sure of his formulz it is now al- 
most necessary that every astronomer should have a development of his own. 
This condition seems to be unfortunate, since it causes a repetition of work that 


should be done once for all. 


U. S. N. Observatory, Nov. 5, 1886. 
Communicated by the Superintendent. 


ON THE FORM AND POSITION OF THE SEA-LEVEL AS DEPENDENT ON 
SUPERFICIAL MASSES SYMMETRICALLY DISPOSED WITH RE- 
SPECT TO A RADIUS OF THE EARTH’S SURFACE. 


By Mr. R. S. Woopwarp, Washington, D. C. 
[CONTINUED FROM VOL. II, PAGE 131.] 

18. We proceed now to the determination of the constants l, of equation 
(6) and U, of (56) and (57). It has already been stated that these constants are 
to be determined from the condition of equality in volumes contained by the dis- 
turbed and undisturbed surfaces, a condition whose analytical statement is 

vsinada =o. 


* Philosophical Transactions, London, 1804. 
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Substituting the vaiue of 7 from (6) in this, there results 


sinada— |, sinada=o, 


whence =} Vsin ada. (58) 
| 


The easiest way to evaluate the integral in this equation is to substitute the 
value of I” from equation (54). We get then at once 


= sin? $3.4 f sin ada 
Oo 


= 4r,hpz sin’ $3, (59) 


since all terms of the series, except the first, vanish in the integration. Ina sim- 
ilar manner it may be shown that 


U, = sin? 47. (60) 


I’, may also be derived thus: For points within the perimeter of the at- 
tracting mass replace |” in (58) by I, of (20), and for points outside that perime- 
ter replace V in (58) by I, of (21). Making these substitutions there results 


* 
V, = 27,hp [fA sin ad a sin ada| (61) 
re) 3 


Substituting the value of /, from (22), 


3 3 
92129 _ cin? 


sin? 7, 1 — da sin’ 7, 


*A more direct process than that followed in the text is indicated thus : — 


7, sinada + | sinada = —]| /, cosa [| + § cosa da, — | /, cosa | + § cosa da. 
° 3 


0 


| 

| 

| 
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AS DEPENDENT ON SUPERFICIAL MASSES, ETC. 13 i 
1 — sin* 7, 
Then the last integral becomes r 
ty ‘ 
4 cos’ $3 sin?7,) (¢2— 1) | 
ty 
2 1 2) (I — 
in which t, = sin 43, 4 
sin 3,3 cos 7, 
— sin’ sin? | 
Substituting these limits in the non-logarithmic part, there results | 


ir 
+ 4 sin? 43 [ cot’ log. 


Integrating by parts all terms of this expression except the first, we get 


lor 
( (1 — sin? 4,3 sin’ 7,)! cot 7; 
sin 43 sin 7, sin d3 
-++ arc sin (sin 43 sin 
4 sin? $3. (sin 3; 


| — cot? 43 arc sin (sin sin J 


This gives I, sin a da = 2 (sin ? — (62) 
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The second integral in (61) becomes, by substitution of the value of /, from (25), 


sin ada 
B 
. d (sin? 4a) 
J (sin? $a — sin? 43 sin’7,) (1 — sin? 47 sin*y,) 


in 4,3 cos 7 sin sin’ 7, cos 7 
hr 


sin 27, — 2 arc sin (sin sin | 


(1 —sin’ $3 sin’ 7,)! sin 7, 
arc sin (sin $3 sin 72) 
sin? 
whence | /, sin ada = 2 (x sin? $3 — sin ,3 + 3 cos 3). (63) 


The sum of (62) and (63) is 2z sin? 4,3, which, substituted in (61) gives the same 
value for I, as (59). 


19. By reference now to equations (3), (6), (20) to (23), and (59), we find for 
the equation of the disturbed surface, when the effect of the rearranged water is 


neglected, 


v= 3 — sin? 13 (64) 


The corresponding expression in polar harmonics is, see equations (56), (57), 
and (60), 


v= fh! f,(cosa) 


Under the assumption that the water covers the whole sphere and is free to 
adjust itself as stated in §17, the equation to the disturbed surface is 


| 

| 

| 

| | 
| 

| 

| 

| 

| 
| 7 
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20. The position of any point on the disturbed surface is thus defined by 
the co-ordinates v and /, or rather v and a, since / is a function of a; v being 
the elevation or depression of the point relative to the spherical surface, and a 
the angular distance of the point from the axis of the attracting mass. / is to 
be computed from (22) or (23), or their equivalents (24) and (25), according as the 
point is within or without the perimeter of the disturbing mass. The functions 
which enter the last term of (66).are given by (45) and (47) respectively. 

The general nature of the disturbed surface, when the effect of the rearranged 
water is neglected, is evident from (64). It is symmetrical with respect to the 
axis of the attracting mass. It lies without or within the spherical surface of 
reference according as / is greater or less than z sin’ 4,3. The values of / fora 
point of the disturbed surface at the axis of the attracting mass, along its border 
and opposite the centre of the mass, are given by equations (26), (28), and (30) 
respectively. If we denote the corresponding values of 7 by the suffixes 1, 2, 3, 
we get 

== 3h (sin — sin? 4,3), 


3 
= 3/ (é — sin? (67) 
a= 
v5 = 3h = (2 sin? 13 — sin? 1,3). 


The meaning of these equations may be most readily understood by refer- 
ence to Figure 3. Thus if the circle E-GH 
represent (in cross-section) the undisturbed 
sea-level surface of the earth, and a stratum 
of matter, as an ice cap, ABDGE be added 
thereto, the new sea-level surface will assume 
the form indicated by the dotted line. The 
values of v, and 7, as shown in the diagram 
are positive, while the value of 7, is nega- 
tive. If on the other hand we suppose the 
space A’S’D'GFE to be filled with matter of 
less density than the average density of the 
earth's crust, as is the case ina lake basin, 
the disturbed or new sea-level surface would 
fall within some portion PFQ of the undis- re 3. 
turbed surface and outside the remaining portion P//Q; i. e. 7, and v, would be 
negative and 7, positive. 


et 


a 
| 


| 
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21. It is of interest to inquire what angular extent of mass will produce nu- 
merical maxima of 7, 7, and 7, supposing the thickness #, and the densities p,, 
and » constant. By means of the usual criteria it is readily found that 

v,—=amaximum for = 60°, 
=a maximum for sin = 2/z, or = 39°32’, (68) 
amaximum for 120°. 

22. A glance at equation (66) suffices to show that the effect of the free wa- 
ter, if it covers the whole earth, is simply to produce an exaggeration of the type 
of surface defined by (64) and (65). The series in the third member of (66) ex- 
pressing this exaggeration is rapidly converging on account of the diminishing 


2 
Pin 
(2i +1) — 3 


which is, is about 4,1, 


‘ 


factor 


222 5 

The essential features of the disturbed surface are, therefore, in any case de- 
fined by (64) or its equivalent (65); and in most cases the effect of the rearranged 
water may be neglected as unimportant or as of no greater magnitude than the 
uncertainties inherent in the data for actual problems. 

23. The equations (67) define the position of the disturbed surface in some 
of its most characteristic points. To define its position at any other point we 
must evaluate the elliptic integral 7, or /, which pertains to such point. These 
integrals have already been expressed, equation (53), in a series of polar harmo- 
nics, which, if more convergent, would suffice for computing /, or 4. It is easy, 
however, to derive more convergent and convenient series than that of (53), and 
this will be the object of the present section. 

First, take 7, of (24). For brevity put 

sin da 
sin 


- 


and 6 = sin 


Then by Maclaurin’s series, or by the binominal theorem, we readily find 


i= 2b — Asin’ 7, — Bsin'y, — Csin’7,—. . .)dy, (69) 


in which A= tu"? (1 — &), 
B= + 20 — 364), 
C= (1 + P+ — 56), 


etc. 


| 

| 

H 

| 

| 

| 

| 
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The even powers of sin 7; may each be expanded in a series of the form, 
in which ¢, d, ¢, etc. are constants. 
But since frcos 2uy dy = O, 
we shall need in these expansions only the values of c. The value of ¢ in the 
expansion of sin’"y is 
2n(2n — 1)(2a#—2). . .(#+1) ( 


Applying this formula, and making the integration in (69), there results 


Hence if we put* 
1+ 36), 
—)( 1+ 564), 


(1 —&)( 5+ oF + 154 + 350°), 
7 + 12h + 186 + 286° + 630), 
Lo = (1 — &) (21 + 358 + 500 + 700 + 1050 + 2310"), 
etc., 
= — gw’ — — gw —. (70) 
This series converges rapidly except for values of w near unity. In an im- 
portant practical application, wherein ;4 = 38°, (70) gives, using 


* The general value of ¢ is 
1.1.3.5... (2% — 3) 


| 
| 


n(n — 1) 
1.1.3.5 (2% — 
2n (2m — 1) (2"—2).... 1) 3°5 ( 7) 163 

n(n — 1) — 2) 


+1.1.3.5...(2#—9)1.3.5 — 


[| —1.3.5... (2% —1) 62 J 


| 

| 

| 

| 
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terms up to that in zw" inclusive, 7, too great by about 5 per cent. for the case 
zw’ = 1. But this is the most unfavorable case, and one, moreover, for which the 


exact value of /, is known from equation (28). 
By a process entirely similar to that followed above, the expansion of (25) 


gives,* writing for brevity, 
v = cosec $a, 


(i+ 1? (1+ iv?) 

+ (5+ 347+ 5°) 

+ (35 + + 200° + 35°) 

+ (63 + + + + 354% + 


If in this expression we put 


= 


k= 4+ t+ + AHH +... 
ky = ts + +.... 

A= te +..., 


we find (71) 
= b=sin}3, a> 


This series converges somewhat more rapidly than (70). For the case in which 
4 = 38°, and for the extreme value w = 1, using terms to that in w~" inclusive, 
(71) gives 7, too small by about 3 per cent. 


* The general term of /, is 
-5...(2%—1) 


+1.3.5...(2"— 3)1.mp? 
n(n—tI 
1.3 ys 
I 2n (2a —1) — 2). . -(#+2) : pan 
+1.3.5...(2%—7)1.3-5 


J 


| 

| 

| 

| 

it 

it: 

} 
245 + 735 

if 

| 

i} 

i 

1) 

if 

| | | 

oes 

} 

| 

|} 
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For points near the border of the disturbing mass, /, may be expressed by 
a more rapidly converging series than (71). Thus from equation (23) 


B 
= __ cos — cos “| 
Sy (: cos — cos a 


Let cos — cos 4 = 2 sin} (a + sin} (a — 3) = 2a. 
Then 
a 
(1 cos —cosa (cosp—cosa)? (cosp—cosa) | 
B 
dp 
‘cos Pp — cos 4 


I sin} (a + 
sin log. sin — 


nef, 
dp 
(cos — cos a) (cos — cos a)*’ 
etc. ; 
1 aX 
J cosp—cosa? sina da’ 
J (cos 
B 
(cos p — cosa) 2sin?al de J 
etc. 


The integrals in the third and higher terms of the above series are thus seen to 
depend on the integral in the second term. Making the requisite differentiations 
we find to terms of the third order inclusive, 


a@cosa — 2sin*a) ] sin }(a + 3) 

8 sin’ a , (72) 


a=sin}(a + sin} (a — 3). 


| 
| 
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24. Having derived the requisite formula for computing the position of any 
point of the disturbed surface, it remains to determine the slope of this surface 
relative to the undisturbed surface. 


Differentiating equation (64) with respect to a, and dividing the result by the 
radius of the undisturbed surface 7,, we get 


3hp dl 


Aa” 


(73) 


This expresses the slope or inclination of the disturbed to the undisturbed sur- 
face ina meridian plane through the centre of the disturbing mass; it also ex- 
presses the deflection of the plumb-line in the same plane. 


In order to apply (73), it is essential to have the general value of 


aa. 
{sin}a)" 
Since (> , 
= w" cot 1a; 
and hence (70) gives 
1g,w') 
+ 2g," 
= — fa + >. (74) 
+ 
Similarly (71) gives 
+ 
al, 
+ 


Equations (74) and (75) will suffice for the computation of d// da except for 
points near to or at the border of the attracting mass. As a approaches equality 
to ,j the above series become less and less convergent, and finally divergent 


| 
| 
| | 
| 
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when «4 = 3, or w= 1. This may be most readily seen by differentiating (22) 
or (23) with respect to 4 and then making 4 = 4._ Thus we find 


3 3 


Likewise the integrals (24) and (25) become after differentiating them with respect 
to « and then making « = j3, 


- 


d/ sec? y, tan 7, d; 


da 1 (sec? $3 + tan? 7,)* 


sec? 72 > 
(sec* $3 + tan* 7) (4) 


f sec’ $3 sec? 7, f sec’ 7, tan’ 7, dy, 


(sec? 43 + tan? 7,)° 1 (sec? $4 + tan? 7,)° 


dl, 


da 


This shows the equality of (4) and (4), since (7) is plainly infinite, its value being 


dr 
sec’ $3 
Oo 
25. This failure of equations (74) and (75) for points at the border of the 
attracting mass arises from the fact that the expressions (20) and (21), though 
very approximate for the magnitude of the potential /, are not sufficiently gen- 
eral to give an accurate value of dV’ / da, or the attraction in the direction of the 
arc 4, for those points. To determine the slope of the disturbed surface at the 
immediate border of the disturbing mass a special investigation is requisite. 


Since by equations (3) and (6) the slope is expressed by 


a I aV 
(76) 


rda’ 


we may derive an expression for the attraction dV /7,da directly. The exact ex- 


t 

du. 
| 


22 WOODWARD. ON THE FORM AND POSITION OF THE SEA-LEVEL 


pression for the horizontal attraction towards the axis of the mass, of any ele- 
ment-mass is, using the same notation as in § 7, 


4r* dr sin® $4 cos? $4 dd cos dh | 
V(r + 477 sin? 
and the integral of this is dV 


Now, as heretofore, let r=nmt+4, 

In addition put y=r—r, 
so that ay = adr; 


y=—?, for r=n, 
=h—v, for r=n+h. 
Also let = sin 34, 


whence d= =r, cos 34 dd, 


Making the substitutions and neglecting terms of the order 


“uv 
and (= ] ’ 
2r, 
the above expression becomes ' 
dy cosh di 
Integrating with respect to 7 and substituting the limits given above, there results 
( (h —wv) dz d= 
If now we suppose the attracted point on the border of the attracting mass, the 
limits of ¢ will be o, and, with sufficient approximation, 27, sin 3 cos A= ¢ cos 4, 
say. Integrating with respect to ¢ and substituting these limits, we get 


cos? ccoss 
p — v) cos 4 log, { + 


2 


] cos 4 di. 


It remains to integrate these last expressions with respect to 4 between the limits 


| 
| 
| 
| 
| 
‘| e 
| 


AS DEPENDENT ON SUPERFICIAL MASSES, ETC. 23 


o and 3x. An application of the formula for integration by parts will readily 
transform them to elliptics, but since their element functions decrease very rap- 
idly from the lower to the upper limit, the following process* will suffice. Con- 


sider the 
23 
fcosi di log, LV («+ + 


Vv 


in which / is such that (. —— }" may be neglected in comparison with unity. 


v 
In the cases we have to consider (44) will not exceed ;}5 if cos 4 = +15, 


or 4 = 89° 25’, about. Then since 


cosh vy 


the above integral becomes 


A A a 


fcos A di log, 2c cos A = log, Adi +- | log, cosh 


= sins [ tog. + log, (1 + sin A) 


+ (sin 4 — 1) log cos 4. 


But since sin 4 is very nearly unity, the last expression reduces to 


The error of this integral arising from the use of 4 instead of }z as the upper 
limit, is less than 


2 
(42 — 4) cos log, LV 4) + 
which, if cos 4 = ;}5 and ¢c cos /4/ v = 10, amounts to about x9\55- 


*Given in a somewhat different form by Helmert in 7heorieen der Hiheren Geodisie, Vol. Il. p. 322. 


| 

i 

| 

log, —I. 
USCS’ 
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For the entire attraction, therefore, of the mass for a point on its border we get 


Finally, restoring in the last expression the value ofc, viz. c = 2”, sin ,3,(76) becomes 


hp 8r,sins h—v 
do ry [tog, h log, 1]. (77) 


26. Thus far the disturbed surface has been referred to a spherical surface 
concentric with the earth’s centre of gravity betore the disturbance arose. In 
determining the effects of the ice-mass in glacial times, this is the proper surface 
of reference, since we wish to know the distortion of the sea-level in those times 
relative to the sea-level in preceding and following epochs. If, however, it is 
desired to consider the joint effect in disturbing the sea-level of existing masses, 
like the continents, on the hypothesis that such masses rest on the surface of a 
centrobaric sphere, a better surface of reference will obviously be the disturbed or 
existing centre of gravity of the earth. The use of the latter centre will require 
a slight modification of the preceding formulz defining the disturbed sea-surface. 

To determine the radial displacement of the earth’s centre of gravity due 
to the addition of such a superficial mass as we have considered, it is only neces- 
sary to equate the statical moment of that mass to the statical moment of the 
earth’s mass, the moment plane being perpendicular to the axis of the disturbing 
mass at the undisturbed centre of gravity of the earth. The moment of an ele- 
mentary ring of angular radius ,’, measured from the axis of the disturbing mass, 
is to our order of approximation, 


hoz sin cos a7’. 
Hence if ¢ denote the displacement sought and J/ the earth’s mass, 
B 
Me = { 2 sin ;3’ cos a3’ 
= sin? 3. 
Therefore, by substitution of the value of “J7 given in equation (3), we find 


= sin? 3. (78) 


Now the elevation of any point of the disturbed surface relative to the sphere in 
the new position will be less than its elevation relative to the sphere in the former 
position by an amount whose value to the proper degree of approximation is 


COS G4, 


¢ 
7 
| 4¢ h—v 
=— 29 h + vlog, 
| 
| 
| 
| 
‘4 
| 
| 
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4a being as heretofore the angular distance of the point from the axis of the dis- 
turbing mass. That is, if 7’ denote what 7 becomes by the change in position of 
the sphere of reference, 


w= 7 —ocosa. 


Hence, by virtue of (64) and (78), we find for the equation of the disturbed sur- 
face when the sphere of reference is concentric with the disturbed centre of grav- 
ity of the earth, 


= 3h ri [2 — sin? 33(1 + cosa cos? ; (79) 


and the slope of the disturbed surface with respect to the surface of reference is 


h + |rsinasin’ 3}. (80) 
27. Thus far the thickness of the disturbing mass has been considered uniform. 

To determine the effect of a mass symmetrical about an axis, but of variable 


thickness, we may proceed thus: firstly, suppose the effect of the re-arranged | 


water neglected. Then the differential of equation (64) with respect to gives 


dv — d(/—z sin* 


This expresses the elevation of the disturbed surface due to an annulus of angu- 
lar radius f, of angular width dj, and of height 4, the density » being uniform. If 
in this equation we make / a function of 3, or write 4 = ¢ (3), and integrate be- 
tween the proper limits, the result will be the elevation of the disturbed surface 
due to a mass whose thickness conforms to the law expressed by ¢ (3). Calling, 
for the sake of distinction, the new value of the elevation of the disturbed sur- 
face z’’, and the proper limits of ,3, 3, and ;%,, the result of this integration is 


Bs 
d(/-—zsin? 43) _.. ,, 
Tim di ¢ (3) (82) 


Vv = 


Secondly, if the effect of the re-arranged water be taken into account, we 
must add to 7” of (82) the following increment obtained from the third member 
of (66):— 


i= 0 I, (cos a) (3) 


i=! 


& 


| 
| 
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The equations (82) and (83) assign the effect of any homogeneous mass 
symmetrical with respect to an axis, subject to the restriction that the maximum 
thickness of the mass may be neglected in comparison with the radius of the 


earth. 
The integral in (82) depends on, and will in general be, no less complex than 


J, which is defined by equations (22) to (25). In the application of (83) it is to 
be observed that by (46) 
dF, (3) 


= f; (cos sin 3. 


SOLUTIONS OF EXERCISES. 


94 


O is the centre of the circumscribed circle of ABC, and D, £, F the middle 
points of its sides. Show that 


OD? + OE? + OF? = 2R! (2R’ —r’), 


where A’, 7’ are the radii of the circumscribed and inscribed circles of the triangle 
of the feet of the altitudes. [R. D. Bohannan.] 
SOLUTION. 


Let // be the orthocentre, and A’, 4’, C’ the feet of the perpendiculars. 
Since OD = cos A, etc., 
+ + OF* = R?(cos’ A + B + cos? C) 
= (1 — 2cos A cos Bos C). 
But R=2R', and wr =A'HcosA 
= BH cos A cos C 
= 2Rcos A cos C; 


OD? + OF? + OF? = = 2R’ — r’). 


{Marcus Baker.] 


ie 
| 
| 
|| 
| 
| 
| 
| 
| 
| 


SOLUTIONS OF EXERCISEs. 27 


97 
In the triangle ABC two lines drawn from C trisect the side A. Given c¢, 
C, and the angle ¢ between the trisecants ; to solve the triangle. 
[Marcus Baker.| 
SOLUTION, 
Let MY, N be the points of trisection. Then the anharmonic ratios of the 
row AMNB and the pencil C-AJ/NB being equal, we have 
sin ACV.sinBCN AM.BN 
sin ACB. sin WON” AB. MN’ 
or, with the notations given, | 
sin x sin vy = }sin Csing, 


where a=ACM, y=BCN. 
We have also at+y=C-—g. 
But since cos (x y) = cos + y) + 2sin rsiny, 


these relations give 
cos — y) = ¢cos(C — ¢) + ¢); 
whence + — 9 is to be found, and thence + and y. A and BP are then found from 
the relations 
_ tan} (x¥—¢) 
tan} (d —NV)= 
tan* 3 (7 + ¢)’ ] 
£. Hendricks. 
98 
THE eccentric anomalies of three points on an ellipse are /,, f., 3. Show 
that the area of their triangie is 


— ps Ps sin hi—Ps. 


J = 2aé sin 
2 2 2 


the centre of its circumscribing circle 


2 
a 2 2 2 
y= sin 231 sin As. 
b 2 2 
and hence show that the centre of curvature of the ellipse at (1, 7) is 


X=>+ > Y= — 
a’ [W. Thornton.]} 


| 

| 
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SOLUTION I. 
The co-ordinates of the three points are 
4, =acosp, 4 = bsinp;; 
Ag =ACOS fo, OSin 
43 =acosp,, 13 = 


— 4, I3 —JI2 
I 43 


cos — cos f, sin fy — sin 
cos cos fp, sin p, — sin 
sin} + p2) sin — f2) cos (fi + pz) sin (fi — 
sin } (pz + fs) sin} fs) cos} ps) sin (p2 — Pps) 
= 2ab sin } — sin } — fy) sin — fr. 


The co-ordinates of the centret are 


= In 


1 
Substituting the values of the co-ordinates of the three points, | 
1 sinf, a? — (a? — &) sin’ 
a=6b 1 sin~p, a*—(a?— sin’ 
I sinf, —(a*— &) sin’ 
I sin’ 
=—dé(a—6)1 sin 
sinp, sin? 
= 6 (a* — &) (sin p, — sin f,) (sin — sin (sin ~, — sin 
= 86 (a* — &) cos} + COS (Pz + fs) Cos + sin — pe) 
sin } — ps) sin — 


*Salmon’s Conic Sections, % 36. 
tSalmon’s Conic Sections, 22 80, *94. 


i 4 

1] 

if ry 

i 
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Similarly 
1 cosp, & + (a?— &) cos’ 
B=a\t cosp, & + (a®*— &)cos*p, 
1 cosp~, &-+ (a? — &) cos’ 
= — 8a (a’ — &) sin} (p, + sin } (P2 + fs) sin (Ps + sin — 
sin — ps) sin } — py). 
Hence, if = a — #, 


+S cos +f) cos} (A +f) +A), 


y=—F sin} + py) sin + sin} (5 + A). 


If the three points coincide, these become 
3 


X= + cos’ p= + 
[Ormond Stone.] 
SOLUTION II. 


24= (45 — (%, — 42) (3 — 
= ab [(sin p, — sin (cos — cos f,) — (sin — sin f,) (cos — cos f,)] 
= 4ab[ — cos }(p, + fp.) sin — sin + sin (Pp, — 
-+ cos + fz) sin — fz) sin + sin} — 


= 4ab sin — ps) sin (P, — fz) [sin } (P2 + fs) cos $ (P, + 
— cos + fs) sin} (p, + 
= 4ab sin } (pz — ps).sin } — fz) sin } (Ps — Py). 
The equations of the perpendiculars to the chords at their middle points are 


»)] 


I —4 + I's) = — [+ —$(42 + 
Eliminating vy, we have 


(42 — 4) (4% + (43 — 42) (43 + 


| 

| 

| 

| 

| 

| 
| 

| 

{ 


30 SOLUTICNS OF EXERCISES. 


whence, after substituting the values of 1, ”, etc., and making suitable reduc- 
tions, we find 


a2 
x= + cos} + ps) cos $(p, + cos (pi + fr). 
In like manner, by eliminating 7, 
sind + sin (Ps + fi) sin} (Pi + pr): 


Suppose (4, 7) and (43, 7) to approach (2, 4), then the limiting position of the 
centre of the circumscribing will be the centre of curvature at (2, 72). But at 
the limit, 


cosp=~, and sinp 
[ Charles Puryear.] 
99 


From a full cask of wine a quantity is taken out at random and the cask 
filled with water, and then a quantity of the mixture is taken out at random and 
the cask again filled with water. What is the probability that the cask now con- 
tains more wine than water ? [Artemas Martin.] 


SOLUTION, 


Let + = quantity of wine taken out at the first drawing, y = quantity of the 
mixture taken out at the second drawing. Then, a being the capacity of the 


cask, a [ -- 4 E —f | == quantity of wine remaining in the cask after the 


a 
second drawing. Putting 

a(1—*] =40, (1) 
I 
we get = 1. 
a 


Now the greater + is, the less the amount of wine left; hence + may have any 


value from 0 to 1. 
I 


is the probability on the supposition that is known. 


y= 


| 
| 
| | 
| 
| 
|. 
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The greatest value y can have is found by putting + = 0 in (1), which gives 
y=ta. Hence the required probability is 


4 
=" [y+ Jatog [1-2] ] 


= 1—4(1 + log 2). [Artemas Martin] 
101 
Express in terms of sin ~' and sin ~'y 
SOLUTION. 


If sin~'x = a and sin~'y = 4, we have 
a+y —tan-1 Sine + sind 
cos a cosé 
= tan~'tan} (a+ 6) =} (sin~'x + sin-'y). 
[ Cooper D. Schmitt.] 


tan 


EXERCISES. 


111 
A circle is inscribed in an isosceles triangle; another in the space between 
the first circle and the vertex, and so on ad infinitum. What is the vertical angle 
when the sum of the areas of the circles is 1 / # of the area of the triangle? 


L. G. Weld. 
112 [ 


To inscribe in the given triangle AAC a triangle POR whose sides QR, RP, 
PQ are perpendicular respectively to BC, CA, AB. [Alfred C. Lane.]} 


fre ja 
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113 
Tue transversal ZN meets the sides AB, AC of the fixed triangle ABC, and 


makes 
AM.AN=BM. CN. 
Find its envelope. [O. Root, Jr.] 
114 
A CIRCLE intersects a conic in four points, P,, 72, P,P, Show that if the 


4-axis be parallel to the axis of the conic the area of their quadrilateral is 


(4 — 4) (1 — J). [R. H. Graves.] 
115 
IN an equilateral triangle inscribe a triangle similar to a given triangle and 
having a maximum area. - [W. M. Thornton.] 


Finp the value of f Qdx where 


Q = cos + cos (a,x + b,)... cos + 
. . and d,,,. . . 6, being constants. [A. Hall.] 


SELECTED. 


117 
Suow that the angle between the tangent to the ellipse @a° + ay’ =a’ 
and the tangent at the corresponding point of the principal circle is given by the 
formula, 


a tan 
(1 —a) + tan*¢’ 
where « = (a — 6)/a is the ellipticity, and ¢ is the eccentric anomaly. 
118 
SHow that the eccentric anomaly which gives the maximum deviation is 
tan~', (1 —a), and find the co-ordinates of the point of contact 7 and the 
equation to the tangent ¢ 


tan 0 = 


119 
Suow that the portion of ¢ intercepted between the axes equals in length the 
sum of the semi-axes a, 4, and is divided at 7 into these two parts. 
120 
Suow that the circumcircle of ¢ry passes through the centre of curvature at 
7, and that the area of the circle of curvature equals that of the ellipse. 
121 
Suow that when the ellipse varies subject to the condition, a + 4 constant, 
the locus of 7 is the hypocycloid, 
ah + yi=(a + OV, 
which touches each ellipse at the point 7 
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